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Abstract. 

We obtain a constraint on the parameters of a generalized cubic Galileon gravity model 
exhibiting the Vainshtein mechanism by using multi-wavelength observations of the Coma 
Cluster. The generalized cubic Galileon model is characterized by three parameters of the 
turning scale associated with the Vainshtein mechanism, and the amplitude of modifying a 
gravitational potential and a lensing potential. X-ray and Sunyaev-Zel’dovich (SZ) obser¬ 
vations of the intra-cluster medium are sensitive to the gravitational potential, while the 
weak-lensing (WL) measurement is specified by the lensing potential. A joint fit of a comple¬ 
mentary multi-wavelength dataset of X-ray, SZ and WL measurements enables us to simulta¬ 
neously constrain these three parameters of the generalized cubic Galileon model for the first 
time. We also find a degeneracy between the cluster mass parameters and the gravitational 
modification parameters, which is influential in the limit of the weak screening of the fifth 
force. 
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1 Introduction 

Modifications of gravity theory is an interesting approach to explaining the accelerated ex¬ 
pansion of the universe. However, any covariant modification of general relativity introduces 
additional degrees of freedom, giving rise to a fifth force. This is strictly constrained by 
gravity tests in the solar system. Solar system experiments [1, 2] are in excellent agreement 
with general relativity, requiring that this additional degree of freedom be hidden on the 
scale of the solar system. Such a process is referred to as a “screening mechanism,” which 
is key for a viable modified gravity model. In general, this screening mechanism works in 
high-density regions where the matter density contrast is nonlinear. However, this does not 
work on large cosmological scales. A screening mechanism that characterizes viable modified 
gravity models is an important feature to be tested with observations. 

The chameleon mechanism [3, 4] is a screening mechanism that works in an f(R) gravity 
model and the chameleon gravity model [5-7]. In these models, a scalar degree of freedom that 
gives rise to the fifth force is screened in a high-density region due to coupling with matter. 
The chameleon gravity model and an f(R) model can be viable owing to the chameleon 
mechanism [8]. The Vainshtein mechanism [9] is another relevant screening mechanism, 
which is exhibited in the Dvali-Gabadaze-Porrati (DGP) model [10, 11], the simplest cubic 
Galileon model ([12-15], see also the appendix), and its generalized version [16, 17]. The 
DGP model is an archetypal modified gravity model developed in the context of the brane- 
world scenario. There are two branches of solutions in the DGP model. The self-acceleration 
branch DGP (sDGP) model [18-20] includes a mechanism to explain self-acceleration in 
the late universe, while the normal branch DGP (nDGP) model [21-23] with a cosmological 
constant is a healthy modified gravity model avoiding the ghost problem [24, 25] . The simplest 
cubic Galileon model is also a typical modified gravity model that explains self-acceleration 
of the universe while avoiding the ghost problem. Our generalized cubic Galileon model is a 
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generalized version of the simplest cubic Galileon model that retains important features and 
contains the DGP models. In these models, a scalar field giving rise to a fifth force is screened 
due to self-interaction on small scales where density perturbations become nonlinear. 

Galaxy clusters provide a unique laboratory for testing modified gravity models ex¬ 
hibiting screening mechanisms, because they are objects on the borderline between linear 
and nonlinear scales, that is, between non-screened and screened scales. The authors of 
[26-28] have investigated a cosmological constraint on the chameleon gravity model using 
galaxy clusters. They put a constraint on the turning scale and the amplitude of the fifth 
force at large scales. The authors of [29, 30] have investigated a constraint on a generalized 
Galileon model exhibiting the Vainshtein mechanism, using an observed weak-lensing profile 
of clusters. They put a constraint on the turning scale and the amplitude of modification of 
the lensing potential. 

The purpose of the present paper is twofold. One purpose is a generalization of the 
methodology for testing a modified gravity model with a galaxy cluster. In the present paper, 
we consider a generalized cubic Galileon model. Within the quasi-static approximation, the 
generalized cubic Galileon model is effectively characterized by three parameters, fj,Q, 
and e. Detailed definitions are given later, but, broadly, fi g and /tl are parameters that 
modify the effective amplitude of the gravitational potential and the lensing potential in the 
non-screened region, while e determines the turning scale from the non-screened region to 
the screened region due to the Vainshtein mechanism. The parameters hq are constrained 
by observations of the gas distribution, in particular an X-ray surface brightness profile and 
the SZ effect [31]. However, the parameter //l is only constrained by observations of lensing 
measurements. Therefore, a combination of observations of gas distribution and the lensing 
signal is essential to put a constraint on the three parameters characterizing the modified 
gravity model. We demonstrate how a combination of multi-wavelength observations of a 
cluster is useful to put a constraint on a generalized Galileon model. 

The other purpose is improvement of the analysis in [27] using new X-ray data [32, 33] 
and lensing [34] observations of the Coma Cluster. In our method of testing gravity with 
a galaxy cluster, the modeling of the gas distribution is important. A basic assumption of 
the model for the gas distribution is hydrostatic equilibrium, that is, a balance between the 
gas pressure gradient force and the gravitational force. In the region where the fifth force is 
influential, the condition of the hydrostatic equilibrium is changed, and the gas density profile 
is modified. However, in general, galaxy clusters are dynamically evolving, and a deviation 
from the equilibrium could be influential. Therefore, we first check the consistency of our 
model by comparing theoretical predictions with various observations of the Coma Cluster, 
including the new X-ray data and lensing measurements. 

This paper is organized as follows: In Section 2, we first review our model for the dark 
matter and gas distribution of a galaxy cluster within Newtonian gravity. We demonstrate 
how well our model fits observations of the Coma Cluster in Section 3. We also validate our 
model against the influence of non-thermal pressure. In Section 4, we introduce a generalized 
cubic Galileon model, and explain the modification of our model by the fifth force of the scalar 
field. In Section 5, we discuss degeneracies on parameters and systematic errors focusing on 
special circumstance of using the Coma Cluster. Section 6 is devoted to a summary and 
conclusions. Throughout this paper, we adopt G m o = 0.3, Da = 0.7 and Hq = 70km/s/Mpc, 
and we follow the convention (—, +, +, +)• 
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2 Modeling of cluster profiles 


We review our model for the dark matter and gas density distribution connecting with obser¬ 
vational quantities, X-ray brightness, the SZ effect temperature profile, and the weak-lensing 
profile (see also [27, 29]). Our model is based on an assumption of hydrostatic equilibrium, 
but we also take the possible influence of non-thermal pressure into account. In this section, 
we first consider the case of Newtonian gravity. 


2.1 Mass profile 

We employ a universal mass density profile, namely the Navarro-Frenk-White (NFW) profile 
[35], motivated by predictions of numerical simulations, 


p(r) 


_ Ps _ 

r/r s (l + r/r s ) 2 ’ 


( 2 . 1 ) 


where p s and r s are the normalization and the scale radius, respectively. The mass within 
the radius r is given by M(< r) = 4-7T j[] drr 2 p(r) = 47 Tp s r^m(r), with rri(r) = ln(l + r/r s ) — 
r/r a /(l+ r/r s ). We introduce the concentration parameter c and the virial mass M v j r instead 
of p s and r s , 


r v i r 



M vir = M(< r vir ) = — rl ir A c p c (z) 


( 2 . 2 ) 

(2.3) 


with the virial radius r v ; r . Here, p c is the critical density at the redshift z, and A c = 100 is 
the overdensity contrast determined by the spherical collapse model [36]. 


2.2 Gas profile with hydrostatic equilibrium 

We first assume hydrostatic equilibrium between the gas pressure gradient and the gravita¬ 
tional force in galaxy clusters as 


1 dPtot _ 
p gas dr dr ’ 


(2.4) 


where p gas is the gas density, Ptot = Pth + Pnth is the sum of the thermal pressure, P t h and 
the non-thermal pressure, P nt h (see below for details), and T is the gravitational potential. 
We employ three assumptions to describe the gas physics. First, the equation of state for 
the gas components, P t h = n gas kT gas , where n gas and T gas are the number density and the 
temperature of the total gas component, respectively. Second, the temperature of electrons 
is the same as that of the gas, T e = T gas . Third, the electron pressure satisfies P e = n e kT gas 
with the electron number density n e = (2 + p)n gas / 5, where p = 0.59 is the mean molecular 
weight. In the present paper, for the electron temperature we assume the functional form 

—&2/61 

, (2.5) 


T e (r) = T 0 


1+[ ~J 
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where To, b±, &2 an d n are parameters. Integrating (2.4) with (2.5), we obtain the electron 
pressure profile 



GM(< r) 
^2 


P e (r) = n 0 T 0 exp 


( 2 . 6 ) 










and the electron number density n e (r) = P e (r)/kT e (r), where no is the normalization pa¬ 
rameter of the electron number density, n e . In deriving equation (2.6), we use the relation, 
Pgas = h m p n gas> where rn v is the proton mass. Equation (2.6) is the case of the absence of 
non-thermal pressure; the case including non-thermal pressure is described below. 

Thus our gas distribution model includes 7 parameters (M v i r ,c,no,7o,6i,&2>n)- Using our 
model of the three-dimensional profiles, we construct the observables for the observations of 
X-ray and the cosmic microwave background (CMB) temperature distortion. The X-ray 
emission from clusters are dominated by the bremsstrahlung and line emission caused by the 
ionized gas. For the X-ray observable, we define the X-ray brightness as Bx = norm/area, 
where norm is the spectrum normalization obtained from XSPEC software [37] adopting 
the APEC emission spectrum [38], and area is the area of the spectrum. The spectrum 
normalization is given by norm oc f n e nudV, where nn = 0.86n e is the hydrogen number 
density and V is the volume of the spectrum. Then, we write the X-ray brightness as 


Bx(r±) 


to -14 

47t(1 + z) 2 


n e (r)nu(r)dz [cm 5 /arcmin 2 ], 


(2.7) 


where r± is the radius perpendicular to the line-of-sight direction, which is related with 
r and z as r = + z 2 . The CMB temperature distortion is caused by CMB photons 

passsing through clusters and being scttered by electrons in clusters, can be expressed as 
the difference between the averaged CMB temperature and the observed CMB temperature, 
ATsz, or ^-parameter, 


y(r_ l) = 


AT S z 

27cmb 


<7t 

m e 


P e {r)dz, 


( 2 . 8 ) 


where Tcmb = 2.725 K is the CMB temperature, or is the Thomson cross section, m e is the 
electron mass. 


2.3 Shear profile by gravitational weak-lensing 

We consider a spatially flat cosmological background, and work with the cosmological New¬ 
tonian gauge, whose line element is written as 

ds 2 = — (1 + 2 T)dt 2 + a(t) 2 ( 1 + 2 T)dx 2 , (2.9) 

where a(t ) is the scale factor, and T and T are the gravitational and curvature potentials, 
respectively. The propagation of light is determined by the lensing potential ( < f> — T)/2, 
which means that the weak-lensing signal is determined by (T — T)/2. For example, the 
convergence is given by 

K - ^ £ d X ^ X ~ X '^ X ' A( 2p )($-T), (2.10) 

where x is the comoving distance and A( 2D i is the comoving two-dimensional Laplacian. For 
the case of general relativity, we may set AT = —AT = 4-/r Ga?p. Then, using the thin lens 
approximation, (2.10) reduces to 

« = (XS ~ XL)XL r d X ' [4v xGp{r')\ al (2.11) 

XS Jo 
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where xl and xs denote the comoving distance between the observer and lens and that 
between the observer and the source, respectively, and a l = 1/(1 + zl) is the scale factor 
specified by the redshift of the lens object z^. For a spherically symmetric cluster, (2.11) is 
represented as 

«(n_) = — / dzp{r) (2.12) 

Jo 


with the physical coordinate r = y + z 2 . We define X c = xs/[4vrG(xs — Xl)xl«l]- We 
then define the reduced shear 


lMrx>= 7+<r±) 


1 - K(r ±)’ 

where 7+(rj_) is the tangential shear, which is related with the convergence as 

7+(ri_) = «(< Cl) ~ «(cl), 


(2.13) 


(2.14) 


with 


2 f r± 

t(< r±) = — / drirj./c(ri). 
r i -to 


For the NFW profile, the convergence is given by [39] as 
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(2.15) 


(2.16) 


(2.17) 


with x = r/r s . 

Here, we assume that the source galaxies have random orientation of ellipticity es, the 
average of which is (es) = 0. When we observe the tangential ellipticity of the source galaxies 
fobs.) the average is given by (e 0 b s .) = 9+- Hereafter, we assume that the redshift of the source 
galaxies is (zs) = 0.6, but the results are not influenced by the redshift of the source galaxies 
for nearby clusters. 


3 Consistency test with Newtonian gravity 

In the present paper, we use Coma Cluster observations. The Coma Cluster is one of the best 
observed nearby clusters, and has redshift z = 0.0236. The X-ray distribution [32, 33, 40- 
47], the SZ effect [48] and the weak-lensing measurement [34, 49] have been reported. These 


- 5 - 




































observations revealed that the Coma Cluster has substructures and orientation dependence 
on the gas temperature profiles. The Coma Cluster is thus an unrelaxed system. However, 
we will show that our model based on the hydrostatic equilibrium fits the data of the X-ray 
brightness profiles [32, 33], the SZ effect profile from the Planck measurement [48], and the 
weak-lensing profile by Subaru observations [34]. In general, the assumption of hydrostatic 
equilibrium holds only at the intermediate region of clusters, because of the cooling of the 
gas at the innermost region and the environmental effects at the outermost region. Then we 
use data points in the range 200 kpc to 1.5 Mpc to get rid of systematic effects from the 
innermost and outermost regions of the cluster. 

In this work, we use the observational data of the XMM-Newton [32, 33], which are 
different from those used in a previous paper [27]. In that paper, the weak-lensing profile is 
not used; only the parameters M v ; r and c are used as a prior profile from [49]. However, use 
of the weak-lensing profile is essential to our analysis of the generalized Galileon model. 

To address the theoretical predictions in the previous section with observations of the 
Coma Cluster, we introduce the chi-squared by summing the chi-squared for each observation 
as 

2 _ 2 i 2 , 2 

XXB+SZ+WL — XXB + XSZ + XWL 

where 

2 (Hx(r,)-Hfff ) 2 

Xxb / > a ^obs. ^2 

2 _ \ - {y{r±,i) - Vi hs -) 2 

XSZ / ^ (A?/0-bS')2 

2 (9+(r±,j) ~g+ b f') 2 

X'WL — 2—J (A g ob ?-) 2 ’ 

are the chi-square values for the X-ray brightness, the SZ effect and the weak-lensing, respec¬ 
tively. We note that the covariance of errors is not taken into account in our analysis and 
leave it for future work to study how the observational systematics affect our analysis. 

We perform an MCMC analysis using modified Monte Python code [50] that employs a 
Metropolis-Hastings [51, 52] sampling algorithm. This analysis includes 5 parameters with 
the chi-squared (3.1), X'bx+sz+wl- We require Gelman-Rubin statistics [53] of 1Z— 1 < 0.001 
for each parameter to ensure convergence of our runs. The black dashed curve in each panel 
of figure 1 shows the best-fit profiles for the Newtonian gravity. The minimum value of the 
chi-squared is Xxb+sz+wl/^- 0 ^- = 58/44, and the 2-dimensional marginalized contours of 
the different combinations between the model parameters are shown in figure 4. 


(3.1) 

(3.2) 

(3.3) 

(3.4) 
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Table 1 . Best-fitting parameters and 1-dimensional marginalized constraints ( 68 % CL) to charac¬ 
terize the gas and lensing profiles. To avoid degeneracy between parameters, we fix b\ and 62 simply. 
Our results do not depend on whether these parameters are fixed or not. This table shows the results 
for the Newtonian gravity (second column) and the generalized Galileon model with all modification 
parameters (third column). The minimum chi-squared and the number of degrees of freedom, d.o.f. 
= (number of data points) — (number of model parameters), are listed at the bottom of each column. 


Parameter 

Newtonian gravity 

Modified gravity (full parameters) 

M vir 

l-OSlaoe x 10 15 Mq 

l-04l°;“ x 10 15 Mq 

c 

0 cn+0.23 

O.OJ —0 23 

0 nt+0.21 

n 0 

6.14±°; 2 ® x 10 “ 3 /cm 3 

6.17±g;§f x 10” 3 /cm 3 

T 0 

6.36ig;^ keV 

6.35^? keV 

bi 

2.6 (fixed) 

2.6 (fixed) 

62 

0.5 (fixed) 

0.5 (fixed) 

n 

0.74t°;“ Mpc 

0. 751^07 Mpc 

e' 

- 

0.43 

Pg 

- 

0.24 

Pl 

- 

0.55 

Minimum x 2 /d.o.f. 

58/44 

57/41 
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Figure 1. Best-fit profiles in the frameworks of Newtonian gravity (black dashed curve) and the generalized Galileon model (red solid curve), 
and the observational results. The best-fit parameters are listed in Table 1. Left panel: The X-ray surface brightness from the XMM-Newton 
observations [32, 33]. The errors bars are composed of the Poisson noise and systematic errors that we here assume as 5%. Center panel: The SZ 
temperature profile from the Planck measurements [48]. Right panel: The weak-lensing profile from the Subaru observations [49]. 





























Non-thermal pressure possibly caused by turbulent gas and bulk motion causes a sys¬ 
tematic error when comparing observations of clusters. The authors estimated the fraction of 
non-thermal pressure in the Coma Cluster [54] , which can be larger than that of the thermal 
pressure by 10 percent. Here, we estimate how non-thermal pressure affects our fitting based 
on an estimation with a numerical simulation. To this end, we estimate the hydrostatic 
masses by comparison with the X-ray brightness and SZ effect profiles of the Coma Cluster. 
Here we define the non-thermal fraction f nti h by / nt h = P n th/(P n th + Pth), where P nt h and P t h 
are the non-thermal pressure and the thermal pressure, respectively. In the case including 
non-thermal pressure, the thermal pressure is replaced by P t h = (1 — ./ n th)Ptot- We consider 
the following non-thermal pressure fraction as a function of the radius, 


/nth (j") 


a nt (l + z)^ 


r 500 


Tint 


M200 


3 x 10 14 M q 


n M 


(3.5) 


which is a theoretical prediction with numerical simulations in Ref. [55, 56]. 00 and M200 

mean the radius and mass at the radius where the matter density in the galaxy cluster 
is 500 and 200 times of the critical density, respectively. In the present paper we adopt 
(a n t, /3 n t, n n t,nM) = (0.18,0.5,0.8.0.2), which are the best-fit values in [56] consistent with 
those in [54]. 

The best-fit profile in the presence of non-thermal pressure is not significantly altered, 
compared with the best-fit profile in the absence of non-thermal pressure. Figure 2 shows the 
enclosed mass profiles as a function of radius. The gray-hatched region is the la uncertainty 
interval for the lensing mass. The blue-solid and red-solid regions show the la uncertainty 
intervals for hydrostatic masses fitted without and with non-thermal pressure, respectively. 
The hydrostatic mass estimates are in good agreement with the lensing mass, regardless of 
the inclusion of the non-thermal pressure components. This shows that our fitting method 
is not affected by non-thermal pressure, so we do not consider the non-thermal effect when 
putting a constraint on the modified gravity in the next section. 


4 Testing the generalized Galileon gravity model 


We here consider the generalized cubic Galileon model, with action given by [57], 


S = J d*xy/=g[G 4 (</>)R + K(cj), X) - G 3 (</>, X)U<!> + 




(4.1) 


where K(cj>,X), G^{(f),X) and G^tp) are arbitrary functions depending on the scalar field 
(f> and its kinetic term X = — (<9</>) 2 /2 and £ m is the matter Lagrangian. This model is 
a non-minimal coupling version of the kinetic gravity braiding mode [58], and a subclass 
of the most general second-order scalar-tensor theory [59-61] with G^x = G 5 = 0 , where 
G 4 X = dG^/dX. The simplest cubic Galileon model is the case with K = —X, G 3 oc X , 
and G 4 = Mpj2, where Mpj = l/( 8 irG) is the Planck mass. The DGP model is originally 
a 5-dimensional brane-world model, however, it can be effectively described as a Galileon 
model. Note that the DGP model has two branches of cosmological solutions, the self- 
accelerating branch (sDGP) model [18, 19] and the normal branch DGP (nDGP) model [21]. 
The relation between the generalized Galileon model and the specific models are summarized 
in the appendix. 

We consider perturbations of space-time metric and scalar field. We choose the New¬ 
tonian gauge for the space-time metric (2.9). Assuming spherical symmetry of the system, 
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Figure 2. Spherical masses enclosed within the radii. The gray-hatched region denotes the ler 
uncertainty interval for the lensing mass determined solely by weak-lensing analysis [49]. The blue- 
solid and red-solid regions denote the ler uncertainty intervals for hydrostatic masses without and 
with the non-thermal pressure component, respectively, determined by our joint-fit method. The 
hydrostatic and lensing masses agree with each other, irrespective of the presence of a non-thermal 
pressure component. 


within the sub-horizon scale with the quasi-static approximation keeping with the Vainshtein 
feature, the equations for the gravity and the scalar field lead to [17] 


d T GM(< r) . dQ 

_ GM(< r) dQ 
dr r 2 dr ’ 

dQ__J_(,_ L 8 GA 2 CM(< r) 
dr 4A 2 ^ V r 3 


(4.2) 

(4.3) 

(4.4) 


where Q(x) is perturbation of the scalar field defined by </>(£, x) = 4>{t){ 1 + Q(x)), and 
M(< r) = 47 t fg dr'r' 2 p(r') is the enclosed mass of the halo within the physical radius 
r. Note that the perturbed values \k, 4* and Q in (4.2)~(4.4) are written in the physical 
coordinate. In (4.2)~(4.4), we introduce free model parameters a, £, £ and A 2 , which are 
determined by the arbitrary functions K, G 3 and G 4 . The expressions for a, £, ( and A 2 are 
given in the appendix A. The explicit expressions for the simplest cubic Galileon, the sDGP 
and the nDGP models are also presented there. 

Here, we define the Vainshtein radius ry as 


ry = [ 8 GA 2 CM vir ] 1/3 


8 Ge 2 M v 

H 2- 

n o 


1 1/3 


(4.5) 


10 













where we define e = \J 7/q A 2 £ using the Hubble constant Hq. For r < ry, the scalar field 
can be negligible compared with the Newton potential, so Newtonian gravity is recovered. 
For r>ry the scalar field cannot be negligible, and we have 


dT _ (1 + C(a + Q)GM(< r) 
dr r 2 

d$ _ (1 + CQGM(< r) 

dr r 2 


(4.6) 

(4.7) 


Thus the gravitational and curvature potentials are modified at r ry. These modifications 
affect both the gas and weak-lensing profiles. 

We next construct observational quantities of the gas and weak-lensing profiles consid¬ 
ering the scalar field. Since gas components feel gravitational force through the gravitational 
potential 'k, the X-ray brightness and the SZ profiles are modified through modification of 
\k. On the other hand, the gravitational lensing is characterized by the lensing potential 
(<k — \k)/2, so the modified lensing potential alters the observed lensing profile. We therefore 
introduce the parameters 


/^G = (a + OC, (4-8) 

^L = -(a + 2£)C, (4.9) 

with which we can write dT /dr ~ (1 + po)GM{< r)/r 2 and d(T — <k)/dr/2 ~ (1 + pi)GM{< 
r)/r 2 at r > ry. 

In the generalized Galileon model, with the use of parameters fj,Q, and e our modeling 
for the electron pressure profile (2.6) and the weak-lensing profile (2.12) are modified as 
follows: 


P e (r) = P 0 exp ( J 


dr 


pm p 
kT e (r ) 


GM(< r) fi G 2 


+ Tl H o r ( 1 ~ l/ 1 + 12e 2 —^|m(r 
4e 2 V I/ Pcor 3 


(4.10) 


2 r c 
<r±) = / 

Jo 


dz 


P( r ) ~ 


PLPcO 

2e 2 


^f + Uc 2 t> (r) j 

io{r) - 3p s r 3 m(r)/r 3 
V 7 ! + 12 e 2 p s r 3 m(r) / p c 0 r 3 


(4.11) 


Since the gas pressure tracing the matter density deceases with the cluster-centric radius 
increasing, the pressure gradient is restricted to dP e /dr < 0. This gives the constraints on 
PG- 

Instead of /tq, Pl and e, we introduce p! G = pg/(1 + |//g|)> Pl = Pl/( 1 + |A*l|) 5 an d 
e' = 1 — exp(—e), which span the complete available parameter space of p' G and p' L in the 
interval [—1,1] and e ' in the interval [0,1], respectively. General relativity is recovered when 
Pg = Ml = 0 or e' —>• 1. Using the same method adopted for Newtonian case, we perform an 
MCMC analysis for the modified gravity model including 8 parameters with the chi-squared 
Xbx+sz+wl j defined by (3.1). 

Figure 5 shows the 2-dimensional marginalized contours of the different combinations 
between the model parameters. The best-fit parameters and their 1-dimensional marginalized 


- 11 - 

















Figure 3. The 68% (dark gray regions) and the 95% CL (light gray regions) 2-dimensional marginal¬ 
ized contours for the generalized Galileon model parameters. The red-plus, green-cross and yellow- 
triangle points present models for the simplest cubic Galileon, the sDGP and the nDGP models, 
respectively. 


68% errors are listed in the Table. 1. The red curve in each panel of figure 1 shows the best-fit 
profile for the generalized Galileon model with the minimum value of the chi-squared/d.o.f., 
^xb+sz+wl/^- 0 -^ = 57/41. These profiles almost overlap with the profiles for Newtonian 
gravity (black dashed curves), which shows that the large deviation from Newtonian gravity 
is rejected. We note that there is no significant difference between the red and black curves 
in the best-fit profiles. There is a slight difference in the shear profiles at the large radius 
r > IMpc, which seems to be originated from the large errorbars of the shear data. 

Figure 3 shows 2-dimensional marginalized contours of the confidence levels for the pa¬ 
rameters n' G , // l and e'. fi' G and are parameters from the modification of the gravitational 
potential and the lensing potential, and e' is a parameter characterizing the Vainshtein radius. 
Large values of /tg and m, are rejected at the 68% confidence level, which indicates that the 
possibility of a large deviation from the Newtonian gravity is ruled out, depending on the 
parameter e. When e is smaller, the Vainshtein radius becomes smaller, we can put a tighter 
constraint on hq and /tl- However, e is large, the Vainshtein radius becomes large, which 
makes difficult to distinguish between the Newtonian gravity model and the modified grav¬ 
ity model due to the Vainshtein mechanism. The red-plus, green-cross and yellow-triangle 
points in figure 3 show the representative models, the simplest cubic Galileon model, the 
sDGP model and nDGP model, respectively, at the redshift z = 0.0236. The parameter 
values for each models are shown in Table 2. 
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Table 2. Values of modified gravity parameters for each model at the redshift z = 0.0236. 


Models 

e(e') 

V G (a*g) 

AT (// L ) 

simplest cubic Galileon 

0.77 (0.44) 

0.77 (0.44) 

1.12 (0.67) 

sDGP 

-0.26 (-0.22) 

0(0) 

0.53 (0.43) 

nDGP 

0.20 (0.18) 

0(0) 

0.10 (0.09) 


In a previous work [29], a constraint only on the parameter space p^, and e is obtained, 
based on the lensing observations. As a other recent related work, Barreira et al. investi¬ 
gated cluster masses and the concentration parameters in modified gravity models from shear 
profiles [62]. They focused their investigation on the mass-concentration relation of 19 X-ray 
selected clusters from the CLASH survey in the simplest cubic Galileon and Nonlocal gravity 
models. They found that the mass-concentration relation obtaining from the shear profiles 
for the cubic Galileon model is the same as those for the ACDM model, but no stringent 
constraint on the modified gravity models is obtained. Unfortunately the constraint obtained 
in the present paper is not very stringent too, but one can find the following possibility. We 
emphasize that models with p l = 0 like the sDGP and the nDGP models are indistinguish¬ 
able with Newtonian gravity in the method based on the lensing observations. On the other 
hand, our method of combining the gas and weak-lensing profiles can solve the problem from 
this degeneracy. Future observations would improve the constraint. 


5 Discussion 


5.1 Degeneracies on parameters 

On the MCMC analysis in the previous section, we do not take the range of e, [0,0.1], 
into account because it is hard to converge the MCMC runs because of degeneracy in the 
parameter space. Here, we treat this parameter region for complemental discussion. 

First, taking the limit e —>• 0, which means the fifth force is unscreened everywhere, the 
solutions of the gas pressure (4.10) and the convergence (4.11) are reduced to 


Pe(r) 

n(rj_) 


Pq exp 


/' 

J o 


dr 


(1 + /r)vw 


pm v 
kT c (r) [ 

DO 

dzp(r 


GM(< r) 


(1 + /*g) 


(5.1) 

(5.2) 


Then, the pressure profile and the convergence profile are simply modified by a factor of 
(1 + pq) and (1 + p \,), respectively. In this case, we have P e oc (1 + pQ)M w - ir m(c) /c 3 and n oc 
p s oc (1 + pi,)M v i r c 3 /m(c), then, there are degeneracies between the parameters, M v i r ,c, pq , 
and pi,. Figure 6 compares the results of the MCMC analysis with fixing e = 0 (dark blue 
region (68% CL) and mid blue region (95% CL)) and the results of the Newtonian gravity 
(dark gray region and mid gray region), which is the same as those of figure 4. The best- 
fit parameters are shown in the Table. 3. The CL contours of the blue regions reflect the 
degeneracy between the parameters, M v j r , c, pq, and pi,. 
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Table 3. Same as Table 1 but for the results of the generalized Galileon model in the unscreened limit 
with only fixing e = 0 (second column), and the case with fixing all the modified gravity parameters 
e' = 0.05, /Xq = 0.2, and ^ = 0 (third column). 


Parameter 

Modified gravity (unscreened) 

Modified gravity (fifth force) 

M vir 

1-26 tlil x !0 15 M e 

0.86±g;gf x 10 15 Mq 

c 

o 70 +O.I 3 

0.55 

q c^+0.24 

n 0 

6T5 toil x 10 -3 /cm 3 

6.20 toil x 10 _3 /cm 3 

T 0 

6.35ig+ keV 

6.36l[]+ keV 

bi 

2.6 (fixed) 

2.6 (fixed) 

b 2 

0.5 (fixed) 

0.5 (fixed) 

n 

0.75ig$ Mpc 

0-75i°; 3 6 M P C 

e' 

0 (fixed) 

0.05 (fixed) 

l l G 

-0.10 

0.2 (fixed) 


-0.05 

0 (fixed) 

Minimum y 2 /d.o.f. 

57/42 

60/44 


Next, we show how the presence of the fifth force affects the parameter estimation. 
For example, the blue confidence contours in figure 7 shows the 68% and 95% confidence 
contours of the case with fixing e = 0.05, /tq = 0.2 and = 0. M v i r and c are different 
from those of the Newtonian gravity (gray regions), but other parameters, no, To and b\, are 
not changed. The minimum value of the chi-squared/d.o.f. in the presence of the fifth force 
Xxb+sz+wl/ d-°.f• = 60/44, which is almost the same as the Newtonian case, despite the 
different cluster parameter, M v j r ~ 0.9 x 10 15 Mq (see Table. 3). This result exemplifies that 
the presence of the attractive fifth force affects the estimation of the NFW parameters, M v j r 
and c. This is understood as the consequence of the degeneracy between the modification 
parameters hq and /j ,l and M v ; r and c. 

5.2 Systematic errors 

We shall discuss possible systematic errors. In our analysis, we have assumed spherical 
symmetry for the matter distribution and an equilibrium state for the gas component of 
the balance between the pressure gradient and the gravitational force and the fifth force 
in the case of its presence. We have demonstrated that non-thermal pressure at the level 
suggested by numerical simulations does not alter our results. A future X-ray satellite, 
ASTRO-H [63], will observe turbulent gas motion in the Coma Cluster in more detail, which 
will be informative regarding our result. However, observations of the Coma Cluster suggest 
substructures [34, 49, 64-67] and orientation dependence [40, 41, 47], so the Coma Cluster 
is not thought to be a relaxed system. Dynamical states of the Coma Cluster would give a 
systematic difference between our results and temperature measurement. Our fitting results 
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show that the temperature of the Coma Cluster is around 6.4 keV (see Table. 1), but this 
result seems lower than those of X-ray observations [32, 33, 40, 41, 45-47], which estimate that 
the temperature of the Coma Cluster is around 8-9 keV. Comparing the mass-temperature 
scaling relation for a sample of relaxed clusters [68] with an X-ray temperature observation of 
the Coma Cluster [32, 33], the observed temperature is higher that the temperature expected 
by the mass. The enhancement is at 3 a level of intrinsic scatter [68]. Similar results of high 
temperatures have also been reported by a comparison with other clusters [69]. Depending on 
the orientation and excluding the central region, the temperature of the Coma Cluster could 
be around 6-7 keV [32, 40], but it is difficult to take this dependence into account. Therefore, 
systematic error of temperature in the Coma Cluster would cause a substantial influence to 
the proposed fitting method. In order to reduce a possible dependence of cluster-dynamical 
states and halo triaxiality, it is of vital importance to increase the number of sampling 
clusters. Ongoing and future multi-wavelength surveys such as the Hyper Suprime-Cam 
(HSC) optical survey 1 , the Dark Energy Survey (DES) [70], the eROSITA X-ray survey [71], 
and the ACT-Pol [72] and SPT surveys [73] will be powerful aids to better constraining the 
gravity model. 

6 Summary and Conclusion 

In this paper, we obtained a constraint on the generalized Galileon model through the Coma 
Cluster observations of X-ray brightness, the SZ effect and weak lensing. We have constructed 
a simple analytic model of the gas distribution profiles and the weak-lensing profile (c.f. 
[27-29]). The fifth force affects not only the gas distribution but also the weak-lensing 
profile. In general, the effects depend on different parameters characterizing the generalized 
Galileon model. These features can be investigated by combination of the observations 
of a galaxy cluster reflecting the gas density profile and the lensing signals. Their multi¬ 
wavelength observations are complementary to each other, and are useful to put a constraint 
on the modified gravity model by breaking the degeneracy between the model parameters. 
Systematic study compiling multi-wavelength datasets for a large number of clusters enables 
us to well reduce the systematic errors and constraints on the modified gravity models. 
However, the degeneracy between the parameters, M v j r , c, hq and /tl, persists in the limit of 
the weak screening of the fifth force, which affects the estimation of the cluster parameters. 
Future and ongoing surveys and their joint analysis would be a powerful aid to obtaining a 
more stringent constraint on modified gravity models. 
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A Definitions of the coefficients 

In this appendix, we summarize the coefficients between the generalized Galileon model and 

the specific models used in section 4 (see also [29, 74]). The coefficients in the perturbation 

equations (4.2)-(4.4) are defined as 


a = ai, 



(A.l) 

(A.2) 

(A.3) 


x2 _ 2/3 0 G^H 


(A.4) 



(A.5) 


where 



(A.6) 


(A.7) 


(A.8) 


(A.9) 


0 — — (f>XG^X + 2 HG 4 + (/>Cr40. 


(A.10) 
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These coefficients are determined by the background solution, which follows: 


2 XK x -K + 6 X<j)HG 3X - 2 XG 3(j> - 6 H 2 G 4 - 6 H^G^ = p m , (A.ll) 

K - 2X (g 3<p + 4>G 3X ) + 2 (3H 2 + 2 H^j G 4 + 2 (<£ + 2H<j?j G H + 4 XG 4H = 0, (A.12) 

where p m is the non-relativistic matter energy density and H = a/a is the Hubble parameter. 
The background equation for the scalar field is written as 

j + 3HJ — P = 0, (A.13) 


with 


J = 4>K X + 6 HXG 3X — 2 §G 3 $, (A. 14) 

P = K (t> -2X (g 3H + + 6 (2 H 2 + p) GV (A.15) 


The simplest cubic Galileon model is defined by G 4 = Mpj/2, K = —X and G 3 = 
(r 2 /Mpi)X , which corresponds to taking c\ = —1 in Ref. [75], and thus the coefficients in the 
perturbation equations are 


a = 0 , 

£ = ^G 3 G 3X <j) 2 (j), 
t = g 3X 4> 2 
k M ’ 

\ 2 _ G 3 X (f> 

P ’ 

P = —l + 2G 3X (ci) + 2 Htfi) — 4:TrG 3 G 3X 0 4 . 
When we adopt the late time de Sitter attractor solution [76], 


H(a) 

Ho 



1 

2 



+ 



2 

+ 4(1 — H m o) 


The combinations ££ and A 2 £ 


are given by 


£C 

A 2 C 


(i-n m )(2-n m ) 
H m (5 - n m ) ’ 
( 2 -n m y 


where f2 m (a) = p m (a)/3Mp l H 2 (a ) 


is the matter density parameter. 


(A.16) 
(A. 17) 

(A.18) 

(A.19) 
(A.20) 


(A.21) 


(A.22) 
(A.23) 
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Within the sub-horizon approximation, the DGP model [18, 19, 22] can be effectively 
described by the coefficients 


£ 

C 

A 2 

P 


- 1 , 

1 

2 ’ 

2 

_ 3/3’ 

r 2 
_ ' C 

3/3’ 

1 ± 2Hr c 



(A.24) 
(A.25) 

(A.26) 
(A.27) 
(A.28) 


where the sign “=t” in /3 represents the case of the sDGP model with ” sign and the nDGP 
model with “+” sign. For the sDGP model, we adopt the self-accelerating background 
solution, which is specified by the modified Friedmann equation in the sDGP model [21], 

H{a) 1 — D m o /fl m o | (1 — D m o) 2 , . 

“ —2— + + —i—• (A - 29) 

and r c = 1/(1 — D m o)F7o- On the other hand, the nDGP model has no self-accelerating 
solution without introducing the cosmological constant [11, 13]. Here we consider the nDGP 
model with introducing a dynamical dark energy component on the brane, which is tuned 
such that the background evolves as in the lambda cold dark matter model [77]. 
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Figure 4. The 68% (dark gray region) and the 95% CL (mid gray region) 2-dimensional marginalized contours for the 5 model parameters, 
M v ir [10 14 Mq], c, no [10~ 3 cm -3 ], To [keV] and n [Mpc]. The rightmost plots of each row show the 1-dimensional marginalized constraints (solid) 
and likelihood distributions (dotted). 
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Figure 5. Same as figure 4 but when including the modification parameters e', /j,' g and fi' L in the MCMC analysis. 














































Figure 6. The dark blue region (68% CL) and the mid blue region (95% CL) are the results of the MCMC analysis for the modified gravity model 
with fixing e = 0. The dark gray region (68% CL) and the mid gray region (95% CL) are the results for the Newtonian gravity (same as the figure 

4). 
































Figure 7. The dark blue region (68% CL) and the mid blue region (95% CL) are the results of the MCMC analysis for the modified gravity model 
with fixing e' = 0.05, /j,' g = 0.2 and = 0. The dark gray region (68% CL) and the mid gray region (95% CL) are the results for the Newtonian 
gravity (same as the figure 4). 



























